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Phomopsis blight in Lupinus albus is caused by a fungal pathogen, Diaporthe toxica. It can invade all plant 
parts, leading to plant material becoming toxic to grazing animals, and potentially resulting in lupinosis. Iden- 
tifying sources of resistance and breeding for resistance remains the best strategy for controlling Phomopsis 
and reducing lupinosis risks. However, loci associated with resistance to Phomopsis blight have not yet been 
identified. In this study, quantitative trait locus (QTL) analysis identified genomic regions associated with 
resistance to Phomopsis pod blight (PPB) using a linkage map of L. albus constructed previously from an F 8 
recombinant inbred line population derived from a cross between Kiev-Mutant (susceptible to PPB) and 
Villi A (resistant to PPB). Phenotyping was undertaken using a detached pod assay. In total, we identified 
eight QTLs for resistance to PPB on linkage group (LG) 3, LG6, LG10, LG12, LG17 and LG27 from different 
phenotyping environments. However, at least one QTL, QTL-5 on LG10 was consistently detected in both 
phenotyping environments and accounted for up to 28.2% of the total phenotypic variance. The results of this 
study showed that the QTL-2 on LG3 interacts epistatically with QTL-5 and QTL-6, which map on LG10 and 
LG12, respectively. 
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Introduction 

Lupinus albus L. is grown as a protein crop in numerous 
Mediterranean-type climatic regions (Baer 2006, Luckett et 
al. 2008, Noffsinger et al. 2006) and is primarily used as 
stock feed for sheep and cattle (Hill 2006, May et al. 1993). 
It has potential for increased inclusion in the human diet 
(Erbas et al. 2005, Mariotti et al. 2002), although allergies 
have been reported (Quaresma et al. 2007). 

There are several important yield-limiting fungal diseases 
of L. albus, including anthracnose (Colletotrichum lupini, 
Thomas et al. 2008), Pleiochaeta root rot (Pleiochaeta 
setosa, Luckett et al. 2009), and Phomopsis blights caused 
by Diaporthe toxica (Cowley et al. 2010, 2012b). In com- 
mercial crops of L. albus, Phomopsis infection is often not 
detected until animals exhibit symptoms of lupinosis after 
grazing on infected stubble and fallen seed after harvest 
(Cowley et al. 2010). Lupinosis is a degenerative disorder 
that causes acute liver damage, brain damage and death 



Communicated by M. Ishimoto 

Received June 25, 2013. Accepted January 20, 2014. 

*Corresponding author (e-mail: raymond.cowley@pioneer.com) 



(Allen et al. 1979, Luduena et al. 1989). It results from the 
consumption of toxic metabolites produced by the fungus 
(Peterson et al. 1987). Other Phomopsis species that infect 
soybean also produce similar toxins (Balducchi and McGee 
1987). 

In Australia L. albus cultivation has been based on the 
Kiev-Mutant and Ultra varieties since the late 1970s 
(Gladstones 1976). These varieties have had adequate re- 
sistance to D. toxica over a long period (Wood and Allen 
1980, Sweetingham et al. 1998). Epidemics of Phomopsis 
blights on L. albus have been reported in South Africa 
(Jaarsveld and Knox-Davies 1974) and Poland (Kochman 
and Kubicka 1974). This study was prompted by an out- 
break of Phomopsis blight on a crop of Kiev-Mutant in 
southern New South Wales (NSW, Australia) in 2004, which 
resulted in the death of animals that grazed on crop residues 
(Cowley et al. 2010). A virulent isolate of D. toxica (isolate 
DAR80114) capable of invading all plant parts of L. albus 
was collected from the 2004 outbreak and has been used to 
identify L. albus genotypes with genetic resistance in stems, 
leaves and pods (Cowley et al. 2010, 2012b, 2012c). From 
previous work it appears that resistance to D. toxica in 
leaves, stems and pods of L. albus is under independent 
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genetic control (Cowley et al. 2010). Principal component 
analysis showed that resistance to Phomopsis stem and leaf 
blight was related. Resistance to seed infection and 
Phomopsis pod blight was also related, but independently of 
stem and leaf blight. Further work using detached plant part 
assays demonstrated that the genotypes often displayed vary- 
ing resistance to Phomopsis leaf and pod blight in L. albus 
(Cowley et al. 2012b). Shankar et al. (2002a) postulated 
that stem and pod resistance in the L. angustifolius variety 
Tanjil is also controlled by different genes. Considerable 
genetic variation in the resistance to D. toxica exists within 
L. albus germplasm in leaves, stems and pods (Cowley et al. 
2012a, 2012b, 2012c), although knowledge is limited on the 
underlying genetics. 

Ethiopian landrace Villi A is resistant to PPB (Cowley 
et al. 2012b). This landrace is also a source of resistance to 
anthracnose in L. albus (Adhikari et al. 2009, Phan et al. 
2007). A recombinant inbred line (RIL) population exists 
between Villi A (anthracnose resistant) and Kiev-Mutant 
(anthracnose susceptible, Thomas et al. 2008) and has been 
used to develop the first genetic and comparative map of 
white lupin (Phan et al. 2007), although the marker density 
was low. Anthracnose resistance, flowering time and seed 
alkaloid have been mapped using this population (Phan et 
al. 2007). Vipin et al. (2013) developed a DArT array of 
L. albus and increased the marker density in the genetic 
linkage map of the same RIL population derived from Kiev- 
Mutant x Villi A. 

The Kiev-Mutant x Villi A population was phenotyed 
for resistance to Phomopsis pod and leaf blight (Cowley 
et al. 2012b) using detached plant part assays. The popula- 
tion is segregated, with resistance in pods but not leaves. 
The nature of the genetic control of resistance to PPB in 
L. albus is unknown, but is hypothesised to be polygenic 
due to the continuous phenotypic variation that exists in 
structured bi-parental mapping populations assessed for 
resistance (Cowley et al. 2012b). In this study, QTL map- 
ping was undertaken to understand the genetics underlying 
resistance and to locate loci associated with resistance to 
PPB inZ. albus. 

Materials and Methods 

Plant materials 

A subset of an F 8 -derived recombinant inbred line (RIL) 
population, comprising 93 lines derived from a cross be- 
tween Kiev-Mutant (susceptible to D. toxica isolate 
DAR80114) and Villi A (resistant to the same isolate), was 
used for QTL identification for PPB resistance. Seeds of the 
mapping population (Phan et al. 2007) were provided by 
Dr. Huaan Yang (Department of Agriculture and Food, 
Western Australia). 

Experimental design and data analysis using multi-phase 
experiments 

Phenotyping assays in plant pathology using detached 



plant parts are multi-phase experimental processes (Brien 
and Bailey 2006, Smith et al. 2006). This involves growing 
plants in field or controlled-environment trials (Phase 1) and 
then subjecting a sample removed from these plants to dis- 
ease assessment, usually under laboratory conditions (Phase 
2). Each phase may be subject to non-genetic sources of var- 
iation. To be able to separate these sources of variation in both 
phases from genetic sources requires a multi-phase experi- 
ment with an appropriate experimental design and statistical 
analysis (Cowley et al. 2012b). To achieve this, separate 
randomization is required for each phase, with additional 
replication in Phase 2 (Smith et al. 2006). 

In each experiment, the spatial arrangement of lupin 
lines in both phases was optimized using DiGGer design 
software (available from http://www.austat.gen.org/files/ 
software/downloads). All data were analysed using 
ASReml-R (Butler et al. 2009) with factors describing the 
spatial arrangement of both phases in the model (Cowley et 
al. 2012b). 

Inoculum preparation 

The single-spore isolate of D. toxica used in this study 
was isolated from an infected commercial crop of Kiev- 
Mutant grown in Tarcutta, NSW, Australia, in 2004 (Cowley 
et al. 2010). The isolate (accession no. DAR80114) was de- 
posited in the Living Culture Collection in the Department 
of Primary Industries, Orange, NSW, Australia. For this 
study, inoculum of the isolate was prepared as described 
previously by Cowley et al. (2010). A spore concentration 
of 5 x 10 6 conidia/ml was used in all experiments described 
below. 

Phenotyping for resistance to Phomopsis pod blight 

The phenotyping experiments and analysis have been 
detailed in Cowley et al. (2012b) and are briefly described 
below. 

Screen-house experiment: In phase 1 of the first pheno- 
typing experiment (hereafter referred to as experiment 1), 
ten seeds from each line from the KievMutant x P27 1 74 
RIL population were grown in 50 cm rows with three repli- 
cates in a field soil screen-house in Wagga Wagga, NSW, 
Australia (latitude: 35.05°; longitude: 147.35°) in a ran- 
domised complete block. The late-flowering lines (n = 26) 
were sown on 18 July 2008, with seeds of the remaining 
mid- to early flowering lines (n = 71) sown 2 weeks later to 
maximise physiological uniformity among the lines at the 
time of removing pods. 

Glasshouse experiment: In phase 1 of the second set of 
phenotyping experiments (collectively referred to hereafter 
as experiment 2), pods were selected from plants grown in 
sandy-loam soil in 125 mm diameter pots in an evaporative- 
cooled glasshouse in Wagga Wagga. Nine seeds of each line 
were sown per pot and later thinned to three plants per pot. 
The entire population was assessed in three separate experi- 
ments, with two pots of both Kiev-Mutant and Villi A sown 
within each replicate as control lines. Pods were collected at 
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physiological maturity, approximately 24 weeks after sow- 
ing. 

The experimental design and analysis of Phase 2 of the 
detached pod assays are described in Cowley et al. (2012b). 
Briefly, in the first phase of the experiment, plants were 
grown in a randomised replicated array to produce pods (as 
described above). The pods were collected when they 
reached physiological maturity, growth stage 4.4 (Dracup 
and Kirby 1996), and transported to the laboratory. Addi- 
tional pods were collected from randomly selected rows to 
provide the duplication needed for the laboratory phase 
(Smith et al. 2006). The pods were immersed in a D. toxica 
spore suspension (5 x 10 6 spores per ml) before being 
placed in humidified tubs in a culture room at 20°C with 12- 
hour fluorescent light. The incubation period was 7 days and 
then the lids were removed to reduce humidity. The pods 
were scored for disease symptoms 10 days after inoculation 
according to a 0 to 9 scale (where, 0 = no symptoms; 9 = 
total pod rot) as described previously (Cowley et al. 2012b). 
A score <3 was regarded as resistant to PPB. 

The three experiments where Phase 1 occurred in a glass- 
house were analysed using meta-analysis performed in 
ASReml-R. The predicted genotype means were then used 
in the quantitative genetic analysis for QTL detection. 

Quantitative genetic analysis 

An integrated genetic linkage map consisting of 441 
markers (220 AFLP, 105 genie and 136 DArT) on 38 link- 
age groups, with a total length of 2,169 cM (Phan et al. 
2007, Vipin et al. 2013), was utilised to identify loci associ- 
ated with PPB resistance. QTL mapping was first conducted 
on the predicted means for the detached pod assays from 
screen-house-grown plants. The predicted means from 
meta-analysis of the glasshouse-grown plants were then 
analysed separately. We performed both one- and two- 
dimensional interval mapping analyses using the Rqtl pro- 
gram of the R Statistical package with an error probability 
of O.001 (Broman and Sen 2009, Sen and Churchill 2001) 
[http://www.rqtl.org/]. Multiple QTL models were scanned 
for epistatically interacting QTLs. The corrected QTL mod- 
el was optimised by an iterative process. The LOD (loga- 
rithm of odds) scores were calculated by Rqtl to provide a 
measure of the likelihood that the observed data were due to 
linkage, compared to the alternative that they were due to 
chance. 

Results 

Evaluation of Phomopsis pod blight resistance 

The Kiev-Mutant x Villi A population had a continuous 
distribution for resistance to PPB in both phenotyping data- 
sets, suggesting that resistance was quantitatively inherited 
(Fig. 1). In experiment 2 (plants grown in a glasshouse), the 
grand mean for PPB phenotyping was higher than in experi- 
ment 1 (plants grown in screen-house) (4.96 and 4.08, re- 
spectively). The predicted mean for Villi A was also higher 



in experiment 2 than experiment 1, but Villi A was resistant 
in both experiments (2.56 and 1.04, respectively). There 
was no difference in the response of the susceptible parent 
Kiev-Mutant. There was a significant variation in disease 
scores between the experiments, although the disease scores 
showed a significant correlation (r = 0.65, P< 0.001, 
Fig. 2). 

Identification of QTL for Phomopsis pod blight resistance 

Results of the QTL mapping of the two phenotyping data 
sets are shown in Table 1. In experiment 1, five QTL ex- 
plaining 8.3% to 36.9% of the phenotypic variance for re- 
sistance to PPB were identified on linkage groups LG3, 
LG6, LG10, LG12 and LG17 (Fig. 3). Of these QTLs, 
QTL-2 on LG3 had the maximum LOD score (9.3) and ex- 
plained 36.9% variance. P27174 contributed alleles for re- 
sistance at QTL-2 and QTL-4, whereas Kiev-Mutant con- 
tributed alleles for resistance at QTL-5, QTL-6 and QTL-7. 
QTL-2 on LG3 showed epistatic interaction with QTL-5 
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Fig. 1. Frequency distributions for phenotyping response to Phomop- 
sis pod blight in a Lupinus albus population of Fg RILs from Kiev- 
Mutant x P27174. In experiment 1, plants were grown in a field-soil 
screen-house. Experiment 2 shows the overall predicted means from 
meta-analysis of three glasshouse-grown experiments. Phomopsis pod 
blight was assessed using detached pods in both experiments using a 0 
to 9 scale, where 0 = very resistant and 9 = very susceptible. 
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4 6 
Phomopsis pod blight rating 
Experiment 1 

Fig. 2. Relationship between phenotyping experiments assessing 
Phomopsis pod blight severity assessed using a detached pod assay 
(0 to 9 scale) in a Lupinus albus Kiev-Mutant x P27 1 74 population 
as F 8 RILs. Dashed lines show the grand mean for each data set. The 
location of the parental lines is shown with a cross for Kiev-Mutant 
and an open triangle for P27174. The regression line is shown with the 
95% confidence interval depicted in grey. 



(LG10) and QTL-6 (LG12) (Table 1). 

In experiment 2, four QTLs explaining 8.0% to 29.8% 
phenotypic variance were identified on LG10, LG27, and 
two on LG3 (Table 1). Kiev-Mutant contributed alleles for 



resistance in QTL-1, QTL-5 and QTL-8, which were detect- 
ed on LG3, LG10 and LG 27, respectively, whereas P27174 
contributed alleles for resistance on LG3 (QTL-3). The 
QTLs with the highest LOD scores (>6.7) were located on 
linkage groups LG3 (QTL-3) and LG10 (QTL-5) and ac- 
counted for >28.2% variance. In addition, they showed sig- 
nificant epistatic genetic effects (Table 1). QTL-5 on LG10 
was identified in both experiments, was mapped within the 
same genomic region (2.5 cM) and was likely to be the 
same QTL (Fig. 3). QTL-3 was mapped within 15.6 cM of 
QTL-2, and QTL-3 and QTL-2 may be the same gene asso- 
ciated with PPB resistance. 

Discussion 

The two phenotyping datasets presented here were obtained 
from two different environments, and the data from each 
were assessed independently. The overall mean of experi- 
ment 2 was higher than experiment 1 ; however, in both data- 
sets the relative ranking of the lines did not differ greatly 
and the correlation between the environments was signifi- 
cant (Fig. 2). 

Consistent identification of QTLs in similar chromosom- 
al positions from both assessments suggests that they might 
be the same genomic regions controlling resistance to PPB 
with isolate DAR80114. QTL-5 on LG10 showed the larg- 
est effect in both experiments. On LG3, QTL-2 and QTL-3 
detected in experiments 1 and 2 were located at different 
marker intervals (within 15.6 cM). It is possible that both of 
these QTLs map to the same genomic region. High resolu- 
tion mapping of the genomic region is required and has not 
yet been performed. The results of this study show that 
QTL-2 and QTL-5 had the largest effect on PPB resistance 



Table 1. Detected QTLs for resistance to Phomopsis pod blight (caused by Diaporthe toxica) using a Lupinus albus Fg-RIL population derived 
from a cross between Kiev-Mutant and P27174 



QTL 


LG 


Map position 
(cM) 


pa 




LOD 


Explained 
variance'' (% Va) 


Additive 
effect' 


Proximal marker 


Experiment 1 (plants grown in field-soil 


screen-house) 














QTL-2 


LG3 


17.5 


<0.001 




9.3 


36.9 


-2.37 


M75E38A83 


QTL-4 


LG6 


61.9 


0.020 


* 


2.1 


10.0 


-0.39 


Lup337 


QTL-5 


LG10 


30.0 


<0.001 




6.4 


27.1 


0.75 


M61E35A66 


QTL-6 


LG12 


10.0 


0.012 


* 


4.3 


19.4 


1.31 


lPms-524453 


QTL-7 


LG17 


69.0 


0.040 


* 


1.8 


8.3 


0.22 


M65E35C315 


QTL-2:QTL-5 rf 


LG3:LG10 




<0.001 


*** 


5.7 


38.6 


-0.49 




QTL-2: QTL-6 


LG3:LG12 




0.007 


** 


3.8 


24.6 


0.48 




Experiment 2 (plants grown in a glasshouse) 














QTL-1 


LG3 


0.0 


0.007 


# * 


2.5 


11.5 


0.58 


lPms-751131 


QTL-3 


LG3 


33.1 


<0.001 


*** 


7.1 


29.8 


-0.37 


M48E38B99 


QTL-5 


LG10 


27.5 


<0.001 




6.7 


28.2 


0.45 


M61E35A66 


QTL-8 


LG27 


12.5 


0.035 


* 


1.7 


8.0 


0.10 


Lup87 


QTL-3 :QTL-5'' 


LG3:LG10 




<0.001 


* * # 


6.0 


25.8 


-0.38 





"*: P < 0.05, **:P< 0.01 and ***: P < 0.001. 

b Proportion of variability explained by the putative QTL. Bold-consistent genomic regions associated with Phomopsis pod blight resistance that 

was scored under glasshouse and field conditions. 
0 Alleles from QTL with positive additive effect value are from Kiev-Mutant. Alleles from QTL with negative additive effect value are from 

P27174. 

d Interaction effect between pairs of QTLs. 



QTLs for Phomopsis pod blight in Lvpinus albus 



87 



o 
m 



cz 
o 

CO 
o 
o 



o 
o 



o 
m 



J] | — | — fi IPms-75113lA 

VlPms-518033 
Ztty-IPms-515JS9 
=H=*-M47E32A266 
M75E38B73 
IPms-533404 
M75E38A83 *■ 
IPms-518139 
O U -lr-+ IPms-527083 
lPms-522689 
IPms-751243 
IPms-751435 
IPms-522741 
IPms-750813 
L M48E38BS9 A 



-Lup215 
-Lup338 



1 



— Lup337A 

— Lup257 
— lPms-517335 



— M77E41B78 
M48E32A71 

•|/^4E32B166 
M75E41B330 

■M76E38A105 

•Lup270 

■M77E32B135 

-M78E32B249 

-6TPBP1 



— M60E38B352 
M60E38A354 
M77E35B177 
M78E32B96 



'v 



-l!e6E32BS9 
~M47E41A1I7 

lPms-515831 
IPms-751978 
IPms-515272 
IPms-516215 
IPms-749857 
•|Pms-51K>D2 
IPms-516618 
IPms-7S0702 
VHPms-751771 
-\\W«1E35A86 «-l 
M59E38A79 
L IPms-524707 
_IPms-749856 
" -IPms-751110 



AnMtS13 

K68E41B72 

IPms-524453 A 





"I" ~ 



-T-IPms-524499 
t-IPms-525247 
flPms-525389 
HPm5-525206 
"lHPms-517115 
V-IPms-632980 
-3flPms-524434 
"l(-[Pms-533193 



-JUNBP 
v Lup125 



— LupS7 A 



- 



IPms-522633 



IPms-522526 

IPnK-627037 

M48E32A188 

M77E32A16S 

CALTL 

LSSRIb 

LSSRla 

DSI 

Lup313 

rWANRI 

^l«5E35C315A 



LG3 LG6 LG10 LG12 LG17 LG27 

Fig. 3. Five linkage groups from the Lupinus albus Kiev-Mutant x P27174 map showing the location of putative QTLs conditioning resistance 
to Phomopsis pod blight caused by Diaporthe toxica. Solid rectangles show QTLs identified from phenotyping experiment 1 where Phase 1 was 
conducted in a screen-house. Grey rectangles show QTLs identified from the glasshouse-grown phenotyping in experiment 2. Triangles indicate 
proximal markers to the putative QTL locations. 



in this population and showed epistatic interaction. 

The positive and negative values representing additive 
effects of the QTLs indicate that both parents contributed 
alleles for resistance to the observed variation in the disease 
scores of the Kiev-Mutant x P27174 RJL population. Kiev- 
Mutant was commercialised in Australia in the 1970s and 
had good field resistance to Phomopsis (Sweetingham et al. 
1998). This resistance was maintained in Kiev-Mutant for 
over 30 years until 2004 when an isolate capable of over- 
coming resistance in Kiev-Mutant was identified (Cowley 
et al. 2010). Consequently, it is not surprising that Kiev- 
Mutant donates alleles for resistance to Phomopsis as this 
variety does have resistance to Phomopsis, but was suscep- 
tible to the isolate in this study. 

None of the QTLs in this study appeared to be co-located 
with loci already published for the same Kiev-Mutant 
x P27174 RIL population for anthracnose resistance (Yang 
et al. 2010) or flowering time (Phan et al. 2007). This obser- 
vation suggests that loci controlling resistance to anthrac- 
nose and PPB are different and not located in a cluster, as 
reported in other key crops such as rapeseed and cereals 
(Delourme et al. 2004, Michelmore and Meyers 1998) 

Anthracnose resistance was identified from Villi A and 
transferred to modern varieties (Adhikari et al. 2009). Two 
QTLs identified on LG4 and LG17 conditioned resistance to 



anthracnose (Phan et al. 2007). Using microsatellite-anchored 
fragment length polymorphism (MFLP), Yang et al. (2010) 
were able to identify three sequence-specific PCR markers 
conditioning resistance to anthracnose, which were desig- 
nated WANR1, WANR2 and WANR3. QTL-7 was detected 
in the vicinity of WANR1 (-13.2 cM), which has been asso- 
ciated with resistance to anthracnose disease in L. albus. In 
the first experiment presented here, the marker WANR1 was 
a flanking marker for QTL-7, located on LG17, at a distance 
of 13.2 cM from the QTL. 

Genetic diversity analysis has shown that P27174 repre- 
sents a 'cluster' of genotypes of Ethiopian origin (Raman 
et al. 2008). It is possible that other Ethiopian accessions 
may have the same loci conferring PPB resistance. Further 
work is required to test linkage disequilibrium for PPB re- 
sistance among L. albus accessions. Further alleles (or loci) 
for Phomopsis resistance in pods, and possibly stems, may 
exist in other landraces from this geographic region. We 
have previously identified resistance to PPB in the Ethiopi- 
an landrace P28507 (Cowley et al. 2012b), suggesting the 
merit of further screening of accessions from east Africa. 
This approach is supported by the numerous Ethiopian ac- 
cessions that have being identified with resistance to an- 
thracnose (Adhikari et al. 2009), although complementation 
studies have not been undertaken. 
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This study established the location of loci and identifica- 
tion of markers closely linked to resistance to PPB in an RIL 
population of L. albus. Validation of the effects of the QTLs 
against different genetic backgrounds is required before 
attempting marker-assisted selection in breeding programs. 
Phenotyping for resistance to PPB is lengthy as resistance 
can only be detected towards the end of a growing cycle 
when large pods are present at the correct physiological 
maturity. This restricts phenotyping events to one (if using 
field-based screening) or possibly two events per year if 
growing plants in controlled conditions per growing cycle. 
The use of molecular techniques to aid breeding is well 
understood and currently employed in lupin breeding for 
anthracnose resistance (Yang et al. 2004, 2010) and 
Phomopsis stem blight resistance in L. angustifolius (Shankar 
et al. 2002b, Yang et al. 2002). Breeding efforts to improve 
PPB resistance in L. albus would be enhanced by marker- 
assisted selection. However, further research is needed to 
convert the AFLP and DArT markers located near the QTLs 
identified in this study into usable PCR-based markers for 
routine marker-assisted selection in L. albus improvement 
programs. 
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